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We systematically investigate the structural, electronic and magnetic properties of Mn„ clusters 
(n = 2—20) within the ab-initio pseudopotential plane wave method using generalized gradient 
approximation for the exchange-correlation energy. A new kind of icosahedral structural growth 
has been predicted in the intermediate size range. Calculated magnetic moments show an excellent 
agreement with the Stern-Gerlach experiment. A transition from ferromagnetic to ferrimagnetic 
Mn— Mn coupling takes place at n — 5 and the ferrimagnetic states continue to be the ground 
states for the entire size range. Possible presence of multiple isomers in the experimental beam has 
been argued. No signature of non-metal to metal transition is observed in this size range and the 
coordination dependence of d— electron localization is discussed. 

PACS numbers: 75.75.+a, 36.40.Cg, 61.46.Bc, 73.22.-f 



I. INTRODUCTION 

The search for magnetic behavior in the transition 
metal clusters is motivated largely by the desire to un- 
derstand how magnetic properties change in the reduced 
dimension. This is a question of considerable technolog- 
ical importance. Several unexpected magnetic orderings 
have already been reported in the reduced dimension. 
This ranges from the prediction of net magnetic moment 
in clusters of nonmagnetic (Rh, Ref. 1) or antiferromag- 
netic (Cr, Ref. 2 and Mn, Refs. 3 and 4) bulk materials 
to the enhancement in magnetic moment in clusters of 
ferromagnetic metals (Fe, Ref. 5 and Co, Ref. 6). 

Manganese clusters are particularly interesting among 
all 3d transition metal elements due to the 4s 2 , 3c? 5 elec- 
tronic configuration in Mn atoms. Because of the filled 4s 
and half-filled 3c? shells and the large energy gap ~ 8 eV 
between these levels and as well as due to the high 4s 2 3d 5 
— * 4s 1 3c? 6 promotion energy of 2.14 eV—, Mn atoms do 
not bind strongly. As a result, Mn2 is a weakly bound 
van der Waals dimer with reported bond dissociation en- 
ergy ranging from 0.1 ± 0.1 to 0.56 ± 0.26 cV depending 
upon the different method of analysis&SiiSiii. This weak 
Mn— Mn bonding has been demonstrated through the 
photodissociation experiments for Mn+ (n < 7) cluster 
cationsiaii^. Consequently, the bulk a-Mn, which has a 
very complex lattice structure with 58 atoms in the unit 
cell, has the lowest binding energy among all the 3c? tran- 
sition metal elements. 

Magnetic properties of manganese clusters are rather 
unusual. According to Hund's rule, the half-filled lo- 
calized 3c? electrons give rise to large atomic magnetic 
moment of 5 [Ib- An early electron spin resonance 
(ESR) study suggested a magnetic moment of 5 /zs/atom 
for very small Mn clusters^. However, Stern- Gcrlach 
(SG) molecular beam experiments on Mns_gg clusters by 



Knickclbcin recently revealed the net magnetic moments 
ranging from 0.4 to 1.7 /is/atom^. This differs both 
from the ferromagnetic (FM) small clusters and from the 
antifcrromagnetic (AFM) bulk a-Mn. This experimental 
results can be explained either way that the individual 
atomic moments are small and ordered ferromagnetically 
or the individual atomic moments remain large but their 
orientation flips from site to site i.e., they are coupled 
ferrimagnctically. In the SG experiment, it is important 
to note the relative decrease in the magnetic moment 
for Mni3 and Mnig, as well as the relatively very large 
experimental uncertainty in the measured magnetic mo- 
ment for MnT^i*. In the present work, we will show that 
the local minima for Mni3 and Mnig arise due to their 
'closed' icosahedral structures, whereas, the large exper- 
imental uncertainty (± 58 % of the measured value) in 
the magnetic moment of Mn7 is plausibly due to the pro- 
duction of different magnetic isomers (in addition with 
statistical fluctuation) in the subsequent measurements. 

Earlier all electron (AE) studies 14 : 15 : 16 found Mn-Mn 
FM ordering for Mn„ (n > 4) clusters, which, in turn, is 
not consistent with the SG experiment. However, Nayak 
et al. first predicted ferrimagnetic ground state for Mni3 
with a total magnetic moment of 33 /xa 1 ^. In consistent 
with the SG experiments, very recent AE studies by Par- 
vanova et nlM^ (n =2-9, 12 and 13) and Jones et almS. 
(n = 5 and 6) reported ferrimagnetic ordering in Mn„ 
clusters. Bricre et alQ used ultra-soft pscudopotentials 
(US-PP) to study the intermediate size Mn„ clusters (n 
= 13, 15, 19 and 23) and found icosahedral structural 
growth with an exception for Mni5. However, their pre- 
dicted magnetic moments differ widely from the exper- 
imental values. This might be attributed to the reason 
that the US-PP may not be appropriate in describing 
the transition metals with large magnetic moments. This 
will be discussed briefly later in the section^J Our main 
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motivation of this work is particularly driven by the SG 
experiments^. Here we shall investigate — (i) The struc- 
tural and magnetic evolution of Mn„ clusters, n=2-20. 
(ii) The sudden drop in the magnetic moment at n=13 
and 19 and the very large experimental uncertainty in the 
measured magnetic moment for M117, and (iii) The possi- 
ble presence of isomers with different magnetic structures 
in the SG experimental molecular beam. 

It has also been found by Parks et al. that the Mn„ 
clusters show a downward discontinuity in their reaction 
rate with molecular hydrogen at n = 16, and this was 
attributed to non-metal to metal transition in Mnig2£. 
But if this is indeed true then there should be a down- 
ward decrease in the ionization potential too. However, 
no such abrupt decrease has been seen in the measured 
ionization potential^. In the present paper, we calculate 
both the spin gaps to investigate this issue. 



II. COMPUTATIONAL DETAILS 

The calculations are performed using density func- 
tional theory (DFT), within the pseudopotential plane 
wave method. We have used projector augmented wave 
(PAW) method^S and Perdew-Bruke-Ernzerhof (PBE) 
exchange-correlation functional^ for spin-polarized gen- 
eralized gradient correction (GGA) as implemented in 
the Vienna ab-initio Simulation Package (VASP)S1. The 
id and 4s electrons are treated as valence electrons and 
the wave functions are expanded in the plane wave basis 
set with the kinetic energy cut-off 337.3 eV. Reciprocal 
space integrations are carried out at the T point. Sym- 
metry unrestricted geometry and spin optimizations are 
performed using conjugate gradient and quasi-Newtonian 
methods until all the force components are less than a 
threshold value 0.005 eV/A. Simple cubic supercells are 
used with the periodic boundary conditions, where two 
neighboring clusters are kept separated by at least 12 A 
vacuum space. For each size, several initial geometrical 
structures have been considered. To get the ground state 
magnetic moment we have explicitly considered all pos- 
sible spin configurations for each geometrical structure. 
For transition metals with large magnetic moments, the 
PAW method seems to be more appropriate (as good as 
the AE calculations) than the US-PP approach^!,. The 
US-PP overestimates the magnetization energies and this 
overestimation is even more large for GGA calculations 
than local spin density approximation (LSDA). This is 
due to the fact that the GGA functionals arc more sen- 
sitive to the shape of the wave functions than the LSDA 
functionals. However, the difference between these two 
methods, US-PP and PAW, are solely related to the 
pseudization of the augmentation charges in the US-PP 
approach, which can be removed by choosing very ac- 
curate pseudized augmentation function, which is then 
computationally expensive. For a better description see 
the Ref. 25 by Kresse and Joubert. 



The binding energy per atom is calculated as, 

E b (Mn n ) = — \E(Mn n ) - n£(Mn)L (1) 
n 

n being the size of the cluster. The local magnetic mo- 
ment A4, at each site can be calculated as, 

M = [ [p T (r) - ^(r)] dr, (2) 
Jo 

where p^ (r) and p- 1 (r) are spin-up and spin-down charge 
densities, respectively and R is the radius of the sphere 
centering the atom. For a particular cluster, R is taken 
such that no two spheres overlap i.e., R is equal to the 
half of the shortest bond length in that cluster. 

III. RESULTS AND DISCUSSIONS 

A. Small Ferromagnetic Clusters (Mri2 - M114) 

The Mn2 dimer is the most controversial among all 
the sizes we have studied. Experiments based on res- 
onance Raman spectroscopy^ and ESR™ observed an 
AFM ground state with a bond length 3.17 A. Ex- 
perimentally the binding energy was estimated to be 
0.44 ± 0.30 eV/atom&aJil. However, previous AE-DFT 
calculationsMiiSiiSiiS, predicted a FM state to be the 
ground state with much smaller bond length (~ 2.60 
A) than that of the experimental value. In agreement 
with these calculations, our present PAW pseudopoten- 
tial calculations with PBE exchange-correlation function- 
als, predict a FM ground state with total spin 10 ps and 
with bond length 2.58 A. Very small binding energy, 0.52 
eV/atom, is essentially the characteristic of a van der 
Waals system. However, the bond dissociation energy 
increases by considerable amount due to the reduction 
of one electron from the Mn2 dimer, i.e. by creating a 
hole in the 4s level. As measured by the photodissoci- 
ation experiment the bond dissociation (Mn + • • • Mn) 
energy increases to 1.39 eV for MnJ cation^S. We find 
that the bond length of Mn2 decreases monotonically as 
the net moment decreases. This is consistent with the 
physical picture that the reduction of the interatomic 
spacing leads to comparatively stronger overlap of the 
atomic orbitals which, in turn, reduces the magnetic mo- 
ment. However, we find Mn2 with total moment 4 and 6 
Pb are not bound. An AF Mn2 is 0.52 eV higher in en- 
ergy with bond length 2.61 A. The present results, along 
with the previous AE-DFT calculationa 14 i 16 i 18 i 19 , do not 
agree with the experimental results. This might be at- 
tributed to the fact that the density functional theory is 
not adequate in treating van der Waals interaction prop- 
erly. However, there is no experimental results available 
in the gas phase Mn2 dimer and the ESR experiment was 
done in a rare gas matrix. Therefore, it is possible that 
the Mn atoms do interact with the matrix, which could 
stretch the bond length and could lead to the observed 
AFM state. 
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FIG. 1: (Color online). Atomic spin ordering of the ground 
state and low-lying isomers for Mn3 and Mn4 clusters. Num- 
bers in the parenthesis represent number of atoms in the clus- 
ter, relative energy to the ground state and total magnetic mo- 
ment, respectively. Bond lengths are given in A. Blue (Gray) 
color represents up or positive and red (dark gray) represents 
down or negative magnetic moment. We will follow the same 
convention throughout. 



The case of Mn.3 is extremely interesting as it could 
have either FM or a frustrated AFM structure. We have 
studied triangular and linear structures. An equilateral 
triangular FM state with total moment 15 hb is found 
to be the ground state with bond lengths 2.74 A and 
binding energy 0.82 eV/atom. The frustrated AFM state 
with total moment 5 [Lb is found to be nearly degenerate 
with the FM ground state (lies only 0.05 eV higher in 
energy). This has an isosceles triangular structure with 
one long and two short bond lengths of 2.50 and 2.45 
A, respectively (Fig^l. The resonance Raman spectra 
studies by Bier et alm^ suggest the ground state to be 
Jhan- Teller distorted D%h structure with an odd integer 
magnetic moment. 

For the Mn4 cluster, we examined three different 
conformations: square, rhombus and tetrahedron. A 
perfect tetrahedral structure with bond lengths 2.7 A 
and binding energy 1.18 eV/atom is the ground state, 
where Mn— Mn coupling is FM with total moment 20 [is 
(FigQJ. Three isomers are found and all of them are 
tetrahedral. A ferrimagnetic state with total moment 
10 hb is only 0.08 eV higher in energy. Another ferri- 
magnetic state with total moment 8 fiB is found to be 
degenerate with the AF state with no net moment and 
they are 0.20 eV higher in energy In all these optimal 
structures the distances between two similar spins (d^\ or 
dn) are larger than those of between two opposite spins 
(dn). Our results are consistent with the previous AE 
calculation a 14 i 16 i 18 i 19 . Ludwig et alM- have studied Mn4 
in solid silicon and observed a 21-line hyperfine pattern 
that not only establishes the four atoms to be equivalent, 
but also the total moment to be 20 [Ib- However, the 
present results can not directly be compared with this ex- 
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FIG. 2: (Color online). Atomic spin ordering of the ground 
and isomeric geometries for n = 5-10. Same ordering has been 
followed as in the Table H| 



periment because of possible Si-Mn interaction and there 
is no available report of magnetic ordering for Mn4 in its 
gas phase. 



B. Mn 5 - Mnio 

As the number of atoms in the cluster (n) increases 
the determination of the structural and magnetic ground 
state becomes a very delicate task as the number of local 
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minima in the corresponding potential energy surface in- 
creases exponentially with n. Therefore, more than one 
geometric and/or magnetic structures of comparable sta- 
bility are possible. In the Fig|2|we depict the atomic and 
magnetic structures for the ground state as well as for 
the closely lying isomers for the size range n = 5-10. As 
it is mentioned earlier, to hit the ground state more re- 
liably, we have studied all possible spin multiplicities for 
several geometric structures for a particular cluster size 
n. Calculated binding energies, relative energies, mag- 
netic moments and two spin gaps are given in the Table 
[|]for the entire size range n = 2-20. 

For the Mn5 cluster, a square pyramid and a trian- 
gular bi-pyramid (TBP) were studied. Transition in the 
magnetic order, from FM to ferrimagnetic, is found. A 
ferrimagnetic TBP is found to be the ground state with 
total spin 3 \ib- The next two isomers are also ferrimag- 
netic in nature with total spins 13 [is and 5 \xb- Both of 
these structures also have TBP structure and lie 0.06 and 
0.07 eV, respectively, higher in energy. The next lowest 
energy arrangement is FM and also has a TBP structure 
with total spin 23 hb and lies 0.19 eV higher in energy. 
Our results are in agreement with the very recent AE 
calculationsi2*i2i2i. However, previously the FM ground 
state was predicted by both Nayak et alM^ and Ped- 
crson et ali&. In the recent SG experiment^, magnetic 
moment was found to be 0.79±0.25 /is/atom, which is 
very close to our predicted value 0.60 ^s/atom for the 
ground state. 

We have investigated both the octahedral and the 
capped trigonal bi-pyramid for Mng cluster. A ferri- 
magnetic octahedral structure with total spin 8 [Lb is 
found to be the ground state with binding energy 1.57 
eV/atom. Another octahedral ferrimagnetic isomer with 
total moment 2 \ib is nearly degenerate (0.02 eV higher 
in energy). The next isomer is also a ferrimagnetic oc- 
tahedra, which possess a total moment of 16 hb and lies 
0.05 eV higher. The next favorable isomer is FM and 
has a total moment 26 [is and is 0.14 eV higher than 
the ground state. In an earlier calculation, Pederson 
et alm& predicted a FM octahedral structure with mo- 
ment 4.33 /xs/atom to be the ground state. However, in 
agreement with the recent AE-DFT calculationo 18 i 19 i 20 , 
present calculation predicts the same ground state and 
isomers. Experimentally measured magnetic moment 
0.55 ± 0.10 /is/atomi lies between that of our predicted 
ground state, 1.33 /zs/atom and the first isomer, 0.33 
Hb /atom, which are almost degenerate. It is possible 
that in the SG experimental beam, multiple isomers were 
produced such that the measured value is almost an av- 
erage of the ground state and the first isomer. 

We have considered pentagonal bi-pyramid (PBP), 
capped octahedron and bi-capped trigonal pyramid as 
the possible candidates for the ground state of Mn 7 . The 
most stable configuration is a PBP structure with fer- 
rimagnetic spin ordering, which has a total moment 5 
Hb- The next two closest isomers also have ferrimag- 
netic arrangements with 7 /is and 3 hb total moments 



and they lie 0.09 eV and 0.20 eV higher than the ground 
state, respectively. Our ground state magnetic moment 
agrees with the earlier calculationsi2ii2iS, though we pre- 
dict isomers with different spin arrangements. However, 
Pederson et al. predicted a FM ground stated. Present 
ground state magnetic moment per atom exactly matches 
with the experimental value, 0.72 ± 0.42 /xs/atomi. We 
would like to note the rather large uncertainty here. We 
argue that the plausible presence of these isomers, with 
total moments 7 [Ib and 3 fiB along with the ground 
state (5 /is), in the SG beam might lead to this high 
uncertainty in the measured value. 

Motivated by our earlier study on Cug2&24 ; we inves- 
tigated three different geometries for Mns, viz, capped 
pentagonal bi-pyramid (CPBP), bi-capped octahedron 
(BCO) and tri-capped trigonal bi-pyramid (TCTBP). 
The BCO structure with total moments 8 /is and 12 
I^b are found to be degenerate ground state. Another 
BCO structure with total moment 10 hb lies only 0.03 eV 
higher in energy. The SG cluster beam experiment has 
reported a magnetic moment of 1.04 ± 0.14 /is/atom 4 , 
which is nearly an average of our predicted values. There- 
fore, our present DFT study together with the experi- 
ment in turn indicate the possible presence of these three 
isomers in the experimental beam with almost equal sta- 
tistical weight. The optimal CPBP and TCTBP struc- 
tures have total moments 14 (is and 12 (is, respectively 
and they lie 0.31 and 0.4 eV higher than the ground state. 
Parvanova et. alJS, found CPBP structure to be the most 
stable, however, their predicted magnetic moment is very 
small, 4 fiB, compared to both of our value and the SG 
experiment. Pederson et. alw& predicted a FM BCO 
structure with moment 32 as the ground state. The 
optimal FM structure for all the three geometrical struc- 
tures have total moment 32 fj,B and lie 1.01, 0.63 and 1.17 
eV higher in energy compare to their respective optimal 
ferrimagnetic structure, respectively for BCO, CPBP and 
TCTBP structures. 

For the Mng cluster, as initial configuration we took 
three stable isomers found for Cug22*24 and a capped and 
a centered antiprism structure. The optimal structure is 
a centered antiprism structure with total moment 7 fig, 
which is in very good agreement with the experimental 
value 1.01 ± 0.10 ^a/atomi. The local magnetic mo- 
ment M. (as calculated form the Eq. 2) shows strong 
environment dependency due to the anisotropy in bond- 
ing. The M. of the highly coordinated central atom is 
very small, —0.22 fis, whereas those of the surface atoms 
are quite high and lie between 3.45 and 3.75 [Lb- Par- 
vanova et ali& have found a similar structure but with 
different spin configuration with total moment 9 /is to 
be the optimal structure. The next two isomers have the 
same geometry and have 13 and 7 hb total magnetic mo- 
ment (Fig(2Jand Table 01. The next isomer is a bi-capped 
pentagonal bi-pyramid, which lies 0.21 eV higher with a 
total moment of 9 (1b- The optimal capped antiprism 
structure lies 0.23 eV higher and has a total moment of 7 
fj, b- Note that, all these structures have 5 spin-up (iVj) 
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atoms and 4 spin-down (N±) atoms. 

Different tricapped pentagonal bi-pyramidal structures 
along with different tetra capped octahedral structures 
were tried as initial structures for Mnxo- Four isomers 
exist with almost the same energy. They lie within ~ 
0.01 eV energy (see Table [Q and Fig. 0). All of these 
have a pentagonal ring and could be derived by removing 
3 atoms from a 13-atom icosahedra. Ground state has a 
total magnetic moment 14 fis, which is very close to the 
SG experimental value, 1.34 ± 0.09 /Lis/atomi 



C. Intermediate size clusters: Mnu - Mri2o 

All the intermediate sized clusters with n = 11-20 
are found to adopt an icosahedral growth pattern. The 
ground state structures and the few isomers along with 
their corresponding spin arrangements are shown in the 
FigCD An icosahedral structure without one apex atom 
is found to be the ground state for the Mn 12 cluster. This 
structure has 7V| = 8 and N± = 4 spin configuration with 
a total moment of 16 [Lb- This value is close to the ex- 
perimentally measured value of 1.72 ± 0.04 /^a/atom^i. 
Recently Parvanova et al. predicted the same geometri- 
cal structure but with comparatively smaller. 1 /is/atom, 
magnetic momenUS. We have found two closely lying 
isomers, which have the same geometrical structure with 
total moments 4 \ib and 18 hb (Table [IJ. Another pos- 
sible icosahedral structure without the central atom lies 
much higher in energy. 

The obvious candidates for the Mni3 cluster are the 
icosahedral, hexagonal close packed (HCP) and cuboc- 
tahcdral structures. The variation of total energy as a 
function of the total magnetic moment is plotted in the 
Fig0]for all these three conformations. The icosahedral 
structure is found to be the ground state with Aj- = 7 spin 
structure. The two pentagonal rings are AFM coupled for 
this structure (Fig|3J). Consequently, the magnetic mo- 
ment is found to be small, 0.23 /^s/atom. This predicted 
magnetic moment is much smaller than those of its neigh- 
boring Mni2 and Mni4 clusters (Fig|SJ), what has been 
indeed predicted by the SG cxperimenli^. Although, the 
present value is much lower than the experimental value 
of 0.54 ± 0.06 ^s/atomi. However, we have found an- 
other icosahedral isomer with magnetic moment exactly 
the same with the experimental value, which lies only 
0.08 cV higher in energy (Table QJ. This structure also 
has Nij = 7. Recently, Parvanova et al. predicted sim- 
ilar magnetic ordering^. The optimal HCP and cuboc- 
tahedral structures (Fig[3J) have relatively higher mag- 
netic moments 9 [Ib (N\ = 7) and 11 ^b (Af = 8) , 
respectively, and they lie much higher in energy, 0.89 eV 
and 1.12 eV, respectively. Nayak et al. first predicted a 
ferrimagnetic state for Mni3. However, their predicted 
magnetic moment is quite high (33 /i_b): all the surface 
atoms are antiferromagnetically coupled with the central 
atomii. 

The ground state of Mni4 is the first complete icosa- 



hedra with a single atom capping. This structure has 
Aj = 9, with a magnetic moment 1.29 ^s/atom. In this 
structure the magnetic coupling between the two pen- 
tagonal rings is FM, which was coupled antiferromagnet- 
ically in the case of Mni3 and consequently, it has small 
magnetic moment. The next two isomers lie very close 
to the ground state: they lie only 0.02 eV and 0.05 eV 
higher and have 1.43 and 1.57 /is/atom magnetic mo- 
ment, respectively. These two isomers (not shown in 
Fig|3J) have the same N\ = 9 spin structure, but with 
their different positional arrangement. The experimen- 
tally predicted magnetic moment, 1.48 ± 0.03 yUs/atom, 
is an average of the ground state and the two isomers 
(Table D), which again indicates that these isomers might 
be produced along with the ground state in the SG ex- 
periment. 

The discrepancy between the present theoretical and 
experimental magnetic moment is rather large for Mh.15. 
The present value is 0.87 ^a/atom, whereas the corre- 
sponding experimental value is 1.66 ± 0.02 /is/atom. 
We have also found several isomers (Table P), but none 
of them are close to the experimental value. The ground 
state and all the closely lying isomers within ~ 0.1 eV 
energy spacing are of derived icosahedral structure. The 
two competing icosahedral structures with 5,1,5,1,3 and 
1,5,1,5,1,2 staking (i.e. without or with the apex atom) 
are possible (Fig[3J). The first kind of structure is found 
to be the ground state, whereas the optimal structure 
for the second kind lies 0.06 eV higher with a magnetic 
moment 13 \xb (Fig|3|). Another structure of the second 
kind is found to be degenerate with this isomer, which 
has a magnetic moment 7 /jb (Table QJ. However, us- 
ing US-PP Briere et al. found a bec structure to be the 
ground state with much smaller magnetic moment, 0.20 
/is/atomSl. In the present case this bec kind of structure 
lies 0.28 eV higher (Table D). 

The same structural trend is observed in the case of 
Mni6, the two different competing geometries have been 
found to be the possible candidates for the Mni6 clus- 
ter. Both of these structures can be derived from the 19- 
atom double icosahedra, which has a 1,5,1,5,1,5,1-atomic 
staking. The ground state has a magnetic moment 1.25 
/is/atom with N-\ = 9 spin structure. This structure has 
5,1,5,1,4-atomic staking: both the apex atoms and one 
atom from the lower pentagonal ring are missing from 
the double icosahedra. The next isomer has the same 
atomic arrangement and is almost degenerate, which lies 
only 0.02 eV higher. This has 1.38 ^/atom magnetic 
moment and the same (iVj = 9) spin ordering. For 
both of these structures the central pentagonal ring is 
antiferromagnetically coupled with the upper and lower 
(incomplete) pentagonal ring. The experimentally pre- 
dicted value, 1.58 ± 0.02 ^Br^ is very close to these 
predicted values and confirms the corresponding ground 
state to be really of this 'strange' staking. The next two 
isomers have a different icosahedral geometry and have 
comparatively smaller magnetic moment, 0.63 (Fig. |3J) 
and 0.50 /xs/atom. They lie 0.06 and 0.1 eV higher, 
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TABLE I: Binding energy, relative energy to the GS (AE — E — Eos), magnetic moment (with a comparison to the SG 
experimenii^-) and different spin gaps, Ai and A2, for Mn„ (n = 2-20) clusters. 



Cluster 



E b 



AE Magnetic Moment Spin Gaps Cluster 





(eV/atom) (eV) 


O-E 


: / atom) 


(eV) 








Theory SG Exppii 


Si 


82 


Mn 2 


0.52 


0.00 


5.00 




0.95 


1.31 




0.26 


0.52 


0.00 




0.47 


0.47 


Mn 3 


0.82 


0.00 


5.00 




0.73 


1.27 




0.81 


0.05 


1.67 




0.63 


0.58 


Mn 4 


1.18 


0.00 


5.00 




0.66 


2.35 




1.16 


0.08 


2.50 




0.45 


0.85 




1.13 


0.20 


0.00 




0.41 


0.41 




1.13 


0.20 


2.00 




1.12 


0.21 


Mn 5 


1.41 


0.00 
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respectively. Both of them have 1,5,1,5,1,3 staking, i.e. 
the 13-atom icosahedra is complete. The two complete 
pentagonal rings are antiferromagnetically coupled. All 
these structures have same number of Nf and N± but 
have two different class of atomic arrangements, which is 
consequently the reason for their large difference in the 
magnetic moment. 

The M1117 cluster follows the same structural trend 
seen in both Mni 5 and Mni 6 . The ground state is a dou- 
ble icosahedra without both the apex atoms, i.e. it has 
5,1,5,1,5 staking. The spin structure is iVj = 10 and the 
central pentagonal ring is AFM coupled with the other 
rings. This structure has a magnetic moment of 1.59 
fiB /atom, which is in excellent agreement with the ex- 
periment, 1.44 ± 0.02 //s/atomi. The next two isomers 
also have the same conformation as well as the same spin 
structure. For this size the structure of the second kind 
i.e. the icosahcdral structure with one apex atom (Fig|3J 
lies rather higher in energy. To our knowledge, there is 
no available report for any other elements where this kind 



of staking has been observed to be the ground state for 
Mnis — Mni7 clusters. 

The Mnig is the 19-atom double icosahedra without 
one apex atom. The predicted magnetic moment is 
1.67 /is/atom, whereas the experimental value is slightly 
smaller, 1.20 ± 0.02 ^/atomi. Next two isomers are 
nearly degenerate and have 1.56 and 1.44 /x^/atom mag- 
netic moment (Table [|J. For all these structures the in- 
tegrated magnetization densities M. for the central pen- 
tagonal bi-pyramid are negative (Fig. 

The double icosahcdral conformation is found to be the 
ground state for Mni9 . The predicted magnetic moment 
is 1.10 /is/atom (Fig. |3J), which is smaller than those 
of its neighboring clusters, what has been predicted by 
the experiments^. Another magnetic structure has been 
found to be degenerate with 1 ^s/atom magnetic mo- 
ment (not shown in Fig. 0) . Both of the structures have 
Ni = 7 and the central pentagonal ring is AFM coupled 
with the other two rings. However, the predicted mag- 
netic moment is larger than the experimentally measured 
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FIG. 3: (Color online) The ground state and a few higher energy structures for the size range n = 11 — 20. Note the grouping 
of the same kind of spins. 



value, 0.41 ± 0.04 /i^/atomi In our case a magnetic 
structure with a magnetic moment 0.47 /x^/atom (N± = 
9), which is very close to the experimentally measured 
value, lies only 0.08 cV higher in energy (Fig[3J). The op- 
timal FCC structure lies much higher, 1.53 eV, in energy, 
which is shown in Fig. 

Two degenerate ground states have been found with 
1.40 and 1.50 /is/atom magnetic moment for Mn2o clus- 
ter. Both the structures have N± = 7 (Fig. EJ) spin con- 
figuration and the conformation can be seen as a singly 
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FIG. 4: (Color online) Plot of the total energy as a func- 
tion of total magnetic moment S(= JVf — iVj.) for icosahedral, 
hexagonal closed pack and cuboctahedral conformations for 
Mni3 cluster. 



capped 19-atom double icosahedra. The central pentag- 
onal ring is antiferromagnetically coupled with the other 
two rings. The predicted ground state magnetic mo- 
ment is larger than the experimental value, 0.93 ± 0.03 
fiB /atom**. However, a different spin structure (N± = 
8) with magnetic moment 0.80 (j,b /atom, which is close 
to the experimentally predicted value, lies only 0.07 eV 
higher (Fig. EJ). 

In the intermediate size range, the grouping of like spin 
atoms i.e spin segregation occurs (Fig. OJ). For a par- 
ticular sized cluster, we find that the ferromagnctically 
aligned atoms have longer average bond lengths^ than 
those of the antiferromagnetically aligned ones. This is 
because of the Pauli repulsion. 



D. Binding Energies 

The size dependence of the ground state binding en- 
ergy for Mn„ clusters (n =2-20) is shown in Fig. [5j Due 
to the lack of hybridization between the half-filled 3d and 
filled 4s states and due to high 4s 2 3d 5 — * 4s 1 3d 6 promo- 
tion energy, the Mn2 dimcr is a weakly bound dimer, 
which is a characteristic of van der Waals bonding 8 910 . 
As the number of atoms in the cluster increases, the bind- 
ing energy increases monotonically due to the increase in 
the s — d hybridization. However, it remains weak as com- 
pared to the other transition metal clusters in the same 
size range. This weak bonding has been demonstrated 
through the photodissociation experiments for Mn+ (n < 
7) cluster cationsii*^. Recently, we have shownSi that if 
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FIG. 5: (Color online) Plot of binding energy per atom as a 
function of cluster size n for the entire size range 2 < n < 20. 
(a) Plot of the same as a function of 1/n for the ferrimagnetic 
clusters, 5 < n < 20 and a linear fit (B.E. = -8.20 ± + 2.80) 
to the data, (b) Plot of second difference, A2E in energy, 
which represents the relative stability. 

an As-atom is doped to the Mn„ clusters, the binding en- 
ergies of the resultant Mn n As clusters increase substan- 
tially due to their hybridized s — d electrons bond with 
the p electrons of As. Similar enhancement in bonding 
has also been seen due to the single nitrogen dopingSi. 

Upon extrapolation of the linear fit to the binding 
energy per atom data to 1/n — > (FigOJa)), we ob- 
tained the binding energy for an infinitely large cluster 
as 2.80 eV, which is very close to the experimental AF 
bulk a-Mn (2.92 eV). It is important here to note the 
kinks observed at n = 7 and 13 in the binding energy 
curve (FigEJ. These kinks represent enhanced stabil- 
ity rendered by their 'closed' geometric structures: Mn7 
is PBP and Mni3 is the first complete icosahedra. If 
this argument is valid then there should also be a kink 
at n = 19, due to the fact it has double icosahcdral 
structure. But we do not see any prominent kink in 
the binding energy curve. So, it will be interesting to 
investigate the second difference in the binding energy, 
A 2 E(n) = E{n + 1) + E(n - 1) - 2E(n), where E(n) 
represents the total energy of an n— atom cluster. As 
A 2 E(n) represents stability of the corresponding cluster 
compared to its neighbors, the effect will be prominent. 
A 2 E is plotted in the Fig^b), where we see a peak for 
Mnig too along with n = 7 and 13. However, in addition, 
without any a priory reason, we observe another peak at 
n = 17, which does not have 'closed' structure (a double 
icosahedra without two apex atoms). 



E. Transition in magnetic ordering 

For very small clusters, n <4, the magnetic coupling 
is found to be FM with magnetic moments 5 fiB /atom, 
which is the Hund's rule value for an isolated Mn atom. 
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FIG. 6: (Color online) Size dependent variation of magnetic 
moment. For the size range 5 < x < 20, it shows excellent 
agreement with the SG experiment. Isomers which lie very 
close to the corresponding GS energy are also shown. 



Although we see that for Mn3 cluster the FM solution 
is nearly degenerate with the frustrated AFM solution. 
The size dependence of the magnetic moment per atom 
is plotted in the Fig|S] We see the transition in the mag- 
netic coupling (from FM to ferrimagnetic) takes place at 
n = 5 and the ferrimagnetic states continue to be the 
ground state for the entire size range n =5-20. FigEl 
shows a very good agreement between experimentally 
measured and our predicted magnetic moments. It was 
seen in the SG experiment that the experimental un- 
certainty in measuring the magnetic moment decreases 
with the cluster size. However, this is not the case for 
M117, for which the measured uncertainty is quite large 
(0.72±0.42 hb, ±58% of the measured value) as com- 
pared to the neighboing sizes. This large uncertainty 
might arise from the presence of isomers with different 
magnetic moments in the SG beam for subsequent mea- 
surements. However, in addition with the statistical fluc- 
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FIG. 7: Plot of Afi Ferli _ F ij) as a function of cluster size n. 
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FIG. 8: (Color online) The s-, p- and d-projected density of states for the central and surface atoms for Mni3 and Mnig in 
their ground state. Gaussian broadening of half- width 0.1 eV has been used. Integrated magnetization density M for each 
atom is given in the box. 



tuation, the above explanation only stands for a plausible 
reason as for all other sizes we did find many isomers with 
different magnetic moments (see Table H] and FigEJ , but 
the corresponding experimental uncertainty is not that 
large. One another striking feature observed in the ex- 
periment is the sudden decrease in the magnetic moment 
at n=13 and 19, compared to their neighbors. Our cal- 
culation reproduces this feature. This is attributed to 
their 'closed' icosahedral structures: first complete icosa- 
hedra for Mni3 and a double icosahedra for Mnig. The 
other geometries studied viz. hexagonal closed packed 
and cuboctahcdral structures for Mni3 and a fee struc- 
ture for Mnig lie much higher in energy, 0.89 eV (9 /j,b), 
1.12 eV (11 Hb) and 1.53 eV (17 hb), respectively, than 
their corresponding ground state. 

In the Fig|5]we have depicted the magnetic moments 
of the very closely lying isomers with their ground state 
(see Table HJ and while comparing those with the ex- 
perimentally observed values, we come to the conclusion 
that for a particular size of cluster, the isomers with 
different magnetic moments are likely to be present in 
the SG cluster beam with a statistical weight and es- 
sentially, the measured moment is the weighted aver- 
age of those isomers. We calculate the energy differ- 
ence between the optimal FM and optimal fcrrimagnctic 
solutions, A^porri-FM) = -E(Fcrri) — E(FM), and plot 
them as a function of cluster size n in the Fig[7] For 
both Mn3 and Mn4 the FM solutions are slightly lower 
in energy than those of their respective optimal ferri- 
magnetic solutions, whereas the optimal FM solutions 
are slightly higher than the corresponding ferrimagnetic 
ground states for Mn5 and Mng. Thereafter, as the clus- 
ter size increases, this energy difference, A_E( Forri _ FM ), 
increases almost monotonically indicating that the opti- 



mal FM solutions become more and more unlikely All 
these optimal FM states have ~ 4 /x^/atom magnetic 
moments. 



F. Coordination and the d-electron localization 

The angular momentum projected local density of 
states (LDOS) show interesting site dependency. The 
s-, p- and d-projected LDOS for the central and surface 
atoms are plotted in the FigJH] for the Mni3 and Mnig 
clusters. We see only d-projected LDOS are significant 
and are of great interest here. The d-projected LDOS 
of both Mni3 and Mnig for the central atoms are broad 
for both majority and minority spin states, which are 
also reflected through their small values of the integrated 
spin densities M. (1.42 and -1.41 /j,b for Mni3 and Mnig, 
respectively). The broadening occurs due to the high 
coordination of the central atom. On the other hand, 
the d-projected LDOS of the surface atoms are rather lo- 
calized and the majority spins are nearly fully occupied, 
which is in agreement with the relatively large local mag- 
netic moments of the surface atoms (3.60 and 3.44 [Lb for 
Mni3 and Mnig , respectively). 



G. Spin gaps: Nonmetal metal transition? 

A spin arrangement in any magnetic clusters is mag- 
netically stable only if both the spin gaps, 
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FIG. 9: (Color online) Plot of spin gaps as a function of 
cluster size n. See Table[I]for the numerical values. 

are positive, i.e. the lowest unoccupied molecular orbital 
(LUMO) of the majority spin lies above the highest occu- 
pied molecular orbital (HOMO) of the minority spin and 
vice versa. We find these two spin gaps to be positive for 
all the clusters (Table |IJ discussed here and are plotted 
in the FigEl <5i an d 82 have local structures, but gener- 
ally decreases slowly as the coordination increases with 
cluster size. Parks et al. found that the Mn„ clusters 
with n < 15 are not reactive towards molecular hydro- 
gen, whereas they form stable hydrides at and above n = 
16, and the reaction rate varies considerably with the 
cluster sizQ 2 ^. They argued it to be attributed from the 
non-metal to metal transition at n = 16. If this is indeed 
the reason, it is likely that the ionization potential would 
show a significant decrease at Mni6, similar to what has 
been observed for free mercury clusters^. Therefore, we 
expect closing up of the spin gaps at n = 16. How- 
ever, Koretsky et al. observed no sudden decrease in the 
measured ionization potential^ and we do not find any 
spin gap closing at Mni6 either. The spin gaps have rea- 
sonable value, Si = 0.33 eV and 82 = 0.22 eV for Mni6 
cluster (FigEJand Table [I}. This abrupt change in the re- 
action rate with H2 at Mni6 is not due to any structural 
change either, as we find all the medium sized clusters 
adopt icosahedral growth pattern and the reason for the 
observed change in the reaction rate remains unknown. 

IV. SUMMARY AND CONCLUSIONS 

We have systematically investigated the structural, 
electronic and magnetic properties of Mn„ (n = 2-20) 



clusters from the first-principles density functional the- 
ory. An extensive search have been made to locate the 
global minima. Due to the intrinsic 4s 2 3d 5 electronic 
structure and high 4s 2 3d 5 — ► 4s 1 3d 6 promotion energy 
Mn-atoms do not bind strongly when they come closer 
to form a cluster. However, binding energy increases with 
the cluster size as the coordination number increases and 
reaches a value 2.37 eV/atom for Mn2o, which is 81 % of 
the bulk value. A magnetic transition from FM to fer- 
rimagnetic ordering takes place at n = 5 and thereafter 
the energy difference between the optimal ferrimagnetic 
and optimal FM structure increases with the cluster size, 
which indicates that the optimal FM states become more 
and more unfavorable with increasing cluster size . How- 
ever, different ferrimagnetic states are possible within a 
small ~ 0.1 eV energy difference and their plausible pres- 
ence in the experimental SG beam along with the ground 
state has been argued. The predicted magnetic moments 
are in agreement with the SG experiment. The sudden 
decrease in the magnetic moment at n = 13 and 19 is 
due to their 'closed' icosahedral structure. It should be 
pointed out here that in the present calculation we as- 
sumed only collincar alignment of spins. However, spin 
canting or noncollinear magnetic ordering is possible in 
small magnetic clusters as it occurs more easily in a low 
symmetry magnetic system^. Icosahedral growth pat- 
tern is observed for the intermediate size range. However, 
to our knowledge, a different kind of icosahedral pack- 
ing have been observed for Mni5 — Mni7 clusters. In 
any particular cluster, the average bond length between 
antifcrromagnctically aligned atoms are 3—8 % shorter 
than that of the ferromagnctically alligned, which can be 
explained in terms of the Pauli repulsion. Spin segre- 
gation has been observed in the intermediate size range. 
The d— electron localization strongly depends on coor- 
dination: localization decreases with the coordination 
number. There is no signature of non-metal to metal 
transition at n =16, which has been predicted^ through 
the downward discontinuity observed in the reaction rate 
with H2. 
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